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Abstract 

In the last years, several models and simulations calculating the radio emis- 
sion from cosmic ray air showers have been developed. However, a number 
of those made conflicting predictions on the pulse shapes and the amplitudes 
of the radio signal. In the scope of this paper, we discuss a detailed compar- 
ison of two independent and complementary theoretical approaches, namely 
MGMR and REAS3. Furthermore, we study the influence of the underly- 
ing air shower models on the predicted pulse shapes and amplitudes and 
show that remaining discrepancies between MGMR and REAS3 are mostly 
determined by the air shower models. With this general agreement, a break- 
through in the understanding of the modelling of radio emission from air 
showers has been achieved. 
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1. Introduction 



Over the past years, the interest in the detection of radio emission from 
extensive air showers (EAS) increased continuously due to its promising po- 
tential and the results derived with experiments such as LOPES jl|, lij and 
CODALEMA js], 0] . To understand the measurements and to study the cos- 
mic rays with the help of radio detection, detailed simulations and a solid 
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theoretical understanding are needed. In the past, the models REAS2 [5|] 
and MGMR ^ made conflicting predictions on the pulse shape and the ra- 
dio pulse amplitudes. While REAS2 predicted unipolar pulses, the pulses 
calculated with MGMR were bipolar and up to a factor of 10 lower (depend- 
ing on observation bandwidth and shower geometry) than the amplitudes 
obtained with REAS2. With REAS3, a missing radiation component re- 
sulting from the variation of the number of charged particles within the air 
shower is taken into account. From that time on, REAS3 predicts bipolar 
pulses 0]. In MGMR, an additional contribution due to the charge excess 
was implemented jsj. These revisions motivate a detailed comparison be- 
tween MGMR and REAS3. 

In the following, an introduction on both models is given to stress that the 
models are completely independent and complementary. 

1.1. A macroscopic approach - MGMR 

When an UHE cosmic ray penetrates the atmosphere, an extensive air 
shower (EAS), a cascade of secondary particles, is created. These particles 
are concentrated in a relatively small plasma cloud that moves with almost 
the light velocity. In this plasma cloud a macroscopic current is induced, 
through the Lorentz force, by the magnetic field of the Earth. The reason is 
that the induced drift velocity, perpendicular to the direction of the initial 
cosmic ray has opposite directions for the electrons and positrons. The knock- 
out of electrons from the air molecules by the high energy photons (through 
Compton scattering) and leptons builds a net excess of electrons that move 
with the same velocity as the plasma cloud. The moving charge excess and 
the induced transverse current can be combined in a four- vector current. 
This four-vector is roughly proportional to the number of charged particles 
in the plasma cloud as the induced drift velocity is varying little with height 
and the electron knockout is a statistical process. As the number of charged 
particles is varying with height the resulting variation of the four-current will 
lead to the emission of electromagnetic radiation. As long as the wavelength 
of the emitted radiation is considerably longer than the size of the charge 
cloud and the (retarded) time-variation of the particle number is larger than 
the oscillation time, all charges contribute coherently resulting in very intense 
radiation. As a result, the frequency spectrum contains information of the 
important length and time scales of the four-current in the EAS j^. The 
polarisation of the emitted radiation offers additional clues on the origin of 
the coherent radiation (sl. [lo|. 



2 



The MGMR-model has its basis in the Lienard-Wiechert potentials from 
classical electrodynamics given by, 



where t' denotes the negative retarded time corresponding to the time the 
signal was emitted from the shower front and V is the retarded distance [sj. 
The volume integral takes care of the particle distributions in the shower 
front. In MGMR the lateral particle distribution is neglected, hence the 
four-current can be written like, 

r{t\x,h) = r{t')5\r)f{h), (2) 

where r = (x, y) is the lateral distance with respect to the shower axis at 
the position z = —ct' + h. Where z is the distance to the Earth's surface 
as measured along the shower path, and h is the longitudinal distance along 
the shower axis with respect to the shower front. From this is follows that 
f{h) describes the normalised longitudinal particle distribution in the shower 
front. In MGMR, this distribution is fitted with a F— probability distribution 
function given by, 

m = h ■ e-^'^/^ ■ -i (3) 

where L denotes the thickness of the pancake which is set to 3.9 m in this 
comparison. The electromagnetic fields are obtained through the standard 
relation, 

for an observer located at x = {x,y,z), and where at t = t' = the shower 
hits the Earth. The velocity of the shower front is c/3 where we set (3 = 1. 

1.2. A microscopic approach - REAS3 

Motivated by the geosynchrotron idea 12j, the C-I-+ based Monte Carlo 
code REASl was published in 2005. REAS calculates the radio emission 
from individual electrons and positrons and superposes their radiation at a 
given observer position. With the updated version, REAS2, it was possi- 
ble to calculate radio emission from air showers which had been simulated 



with CORSIKA 13 . The air shower information derived with CORSIKA 
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is saved in histograms by an interface program COAST. However, it turned 
out that REAS2 (and all other microscopic models at this time) was missing 
a radiation component coming from the variation of the number of charged 
particles. With the latest version REAS3 0], these contributions are now 
taken into account. 

In REAS3, the tracks of the single particles are described by straight seg- 
ments joined by kinks. In each kink, an acceleration process is taking place. 
The radio emission arising from these acceleration processes is calculated us- 



ing the endpoint formalism [IJ]. With this approach all emission related to 
the acceleration of charged particles is taken into account automatically as 
long as the underlying particle movement ist described with sufficient preci- 
sion. Hence, no assumptions on the mechanism of radio emission need to be 
applied. Since also, the radiation at the beginning and the end of the tracks 
is calculated, the radiation due to the variation of number of charges is con- 
sidered as well as the radiation due to the deflection of the particles in the 
Earth's magnetic field. To calculate the emission in each kink, the radiation 
formula derived from the Lienard-Wiechert potentials is used as the starting 
point. With this, the radiation is: 

fx [{f— (5) X (5] 





E(x,t)dt 

' " 'l-/3-f)3/? 

fx {fx 132) 
(1 - /SV) 

where /3i corresponds to the particle velocity before and P2 after the kink 
while /3 = v{t)/c is given by the particle velocity. Furthermore, e indicates 
the particle charge, R{t) = \R(t)\ denotes the vector between particle and 
observer position, f{t) = R{t)/R{t) is the line-of-sight direction between 
particle and observer, and 7 is the Lorentz factor of the particle. The index 
"ret" means that the equation needs to be evaluated in retarded time. 
For a starting-point of the track, only the term with P2 ■ c gives a contribution 
since the velocity before the kink is zero, and vice versa, at the end of the 
track, the contribution is given by the term with /3i ■ c. 

2. Methodology 

To make sure that the results of both models are indeed comparable, a 
set of prototype showers has been defined for the simulations. All of these air 
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Table 1: Overview of the parameters set in MGMR-vl.6 used for the eomparison 
with REAS3. These parameters have been determined from Monte Carlo simulations. 



Parameter (name) 


Value 


Comment 


drift veloeity Vd 


0.025 -c 


c = speed of light 


pancake thickness L 


3.9 m 


assumed as constant 


fraction of charge excess 


0.23 


assumed as constant 



showers have a fixed geometry with a specific energy for the primary particle 
and the observer positions are well-defined. For each prototype shower, a 
set of CORSIKA showers has been simulated. To avoid shower-to-shower 
fiuctuations, one typical shower was selected, i.e. one which has a shower 
maximum Xmax close to the mean Xniax- The chosen hadronic interaction 



models were QGSJetII.03 [15|, ll6| and UrQmdl.3.1 [17|, ll8|. The values for 



the magnetic field have been set to the Argentinian magnetic field at the site 



of the Pierre Auger Observatory [19[, i.e. to a magnetic field of 0.23 Gauss 
with an inclination of —37°, except for the showers considered in section 
13.31 The observer positions were all set to 1400 m above sea level, which 
corresponds to the site of the Pierre Auger Observatory. To calculate the 
radio emission from air showers with MGMR, some parameters need to be 
set which have been determined from Monte Carlo simulations. For this 



comparison they are listed in table [T] and published in 10 



To perform the MGMR simulations, we used the longitudinal profiles of the 
CORSIKA simulated air showers while in REAS3 the histogramed informa- 
tion from the same air shower was taken into account. This ensures that a 
similar parametrisation of the air shower is used for the MGMR simulations 
and thus allows a direct comparison of the results. The difference between 
both, however, is the level of detail in which they take the air shower prop- 
erties into account. Furthermore in MGMR, the lateral distribution of the 
particles in the shower pancake is not taken into account. In section HI the 
differences in the air shower models are discussed in more detail. 
Parametrisations which will not be considered in this article are briefly sum- 
marised in the following: 

In REAS3, the fraction of the charge excess with respect to the total number 
of charges increases with the shower development as given by the longitudi- 
nal development in CORSIKA. In MGMR, the ratio is set to 23% for the 
whole air shower which is comparable to the realistic net charge distribution 
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averaged over the longitudinal development of the air shower. However, this 
leads to differences in the amount of charge excess. 

Furthermore, the velocity of the pancake is set to the speed of light, i.e. 
/3 = 1, in the case of MGMR (cf. section [LT]) . The assumption here is that 
due to the constant creation of high energy electrons and positrons in the 
shower front, the macroscopic particle distribution has the speed of light. A 
direct consequence is that at the back of the pancake a constant flow of low 
energetic particles will be left behind. In REAS3, the particle velocities are 
taken individually according to the histograms. 

An unknown influence might arise due to the underlying atmosphere models 
used by REAS3 and MGMR which are not equal. In REAS, the same atmo- 



sphere models are implemented as given in CORSIKA [13|. For the present 
prototype showers, the US standard atmosphere was used. The atmosphere 
model implemented in MGMR, follows an exponential function for the whole 
air shower development: 



1000.0 g 

X [h\ = ■ exp 

cos a cm"' 



log(0.68) h' 
4000.0 m 



(6) 



where 6 is the zenith angle and h is the height above sea level from where the 
signal is emitted. To calculate the radio emission in REAS3 with a modified 
atmosphere, the atmosphere model would have to be adapted in CORSIKA 
as well, since this influences the air shower development. Thus, the modified 
atmosphere has not been considered for this comparison. 
Note, that for the comparison in this article, the refractive index of the 
atmosphere was set to unity for both models since we want to achieve a 
fundamental understanding of the differences and similarities of MGMR and 
REAS3. Considering the refractive index, the complexity of the radio emis- 
sion mechanism would highly increase. 



3. Direct Comparison 

In table El the set of selected prototype air showers are listed where the 
angles correspond to the coordinate system given by CORSIKA, i.e. an 
azimuth angle of 0° corresponds to an air shower coming from south and 
an azimuth of 90° corresponds to an air coming from east. The set contains 
vertical air showers with various primary energies, one inclined air shower and 
vertical air showers where the magnetic field first was chosen perpendicular 
to the shower axis and second parallel to the shower axis. One shower was 
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Table 2; Overview of the set of prototype air showers whieh were simulated for a 
detailed eomparison between REAS3 and MGMR. An inclination of —37° corresponds 
to the magnetic field geometry of the Pierre Auger Observatory. Note, that the 
azimuth angle is given in the coordinate system of CORSIKA, i.e. the inclined air 
shower is coming from south-east. 



Energy 


Zenith, 
Azimuth 


Strength of B 


Inclination of B 


10^^ eV 


0°, 0° 


0.23 Gauss 


-37° 


10^^ eV 


0°, 0° 


0.23 Gauss 


-37° 


IQi^eV 


0°, 0° 


0.23 Gauss 


-37° 


IQi^eV 


0°, 0° 


0.23 Gauss 


0° (horizontal) 


10" eV 


0°, 0° 


0.23 Gauss 


90° (vertical) 


10" eV 


0°, 0° 


0.0 Gauss 




10" eV 


50°, 45° 


0.23 Gauss 


-37° 



calculated in the absence of any magnetic field to compare the emission of 
non-geomagnetic radiation. This section focuses on three different points: 
vertical air showers (section 13. ip . an inclined air shower (section 13. 2p and 
special magnetic field configurations (section [373|l . 

3.1. Emission from vertical air showers 

To compare the emission from a vertical air shower, three different pri- 
mary energies were simulated. The primary energies were chosen as 10^^ eV, 
10^^ eV and 10^^ eV since these are the typical energies measured by the ra- 



dio enhancement of the Pierre Auger Observatory, AERA (20|. The higher 
the energy of the primary particle Ep, the more electrons and positrons are 
generated in an air shower. A rough estimation for the number of charged 
particles in the shower maximum is given by A'max = Ep/GeV j2l[. The 
scaling of the number of particles with the primary particle energy is directly 
reflected in the scaling of the field strength with particle energy for coherent 



radio emission, i.e. A'max oc E^^^ (cf. section 4.4 in 22|). This increase of 
the amplitudes with larger primary energy is true for the MGMR model as 
well as for the REAS3 simulation as shown in figure [TJ It is obvious that the 
height of the radio signal approximately (but not exactly) scales linearly with 
the energy. This is well understood: the position of the shower maximum 
is deeper in the atmosphere for higher primary particle energies and thus. 



the lateral distribution of the radio signal gets steeper [23| . Combining the 
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Figure 1: Comparison of the east-west polarisation component emitted by vertical air 
showers with three different primary energies at an observer position 200 m north from 
the shower core: 10^^ eV (solid red), 10^^ eV (dashed blue) and 10^^ eV (dotted magenta) 
for REAS3 (left) and MGMR (right). The pulse for 10^^ cV is muhiplied with 0.1 and for 
lO^^eV with 0.01 to allow a better comparison within the same plot. REAS3 and MGMR 
obtain similar results. 



dependence of the primary particle energy and the position of the shower 
maximum, the dependence of the field strength on the primary particle en- 
ergy is still describable by a power-law. 

In fig. [H the pulses of the vertical air showers with different primary 
energies are shown. Since the characteristics of the radio signal are mostly 
unchanged with higher energies 2J] and thus the result of the comparison 
between both models are not influenced by the choice of energy, the following 
comparison is concentrated on the vertical air shower with primary energy 
lO^^eV. 

Figure [1] also shows that the numerical noise level of REAS3 is somewhat 
higher than the noise level of the MGMR simulation. This effect is mostly 
relevant for near-vertical showers as discussed in this section. The figures 
of the inclined air shower illustrate this (cf. section 13. 2p . In the MGMR 
model, the motion of particles is averaged at the beginning of the calculation 
of the coherent emission and the corresponding electric field is calculated at 
the end. Thus, the result of the MGMR model is less affected by numerical 
noise. 

Comparing the raw (unlimited bandwidth) radio pulses from the vertical 
air shower in figure El it is obvious that MGMR as well as REAS3 predict 
bipolar pulses. Furthermore, the amplitudes agree within a factor of ~ 2, 
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Figure 2: Comparison of the west polarisation component emitted by a vertical air shower 
with a primary energy of 10"'^'' eV for REAS3 (left) and MGMR (right). The figures show 
pulses for observers at different lateral distances to the shower core. With increasing 
distance, the results of both models converge. 
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Figure 3: Comparison of the frequency spectra for REAS3 (thick lines) and MGMR (thin 
lines) for a vertical air shower with a primary energy of 10^^ eV. The total spectral field 
strength is shown for observers at different lateral distances from 100 m up to 800 m. 



recalling that the difference between the maximal field strengths of the raw 
pulses of the previous versions of both models was a factor of ten. Only close 
to the shower core, the deviations get larger. This effect can be seen in all 
other figures presented in this section as well. The remarkable agreement for 
larger distances is visible in the frequency spectra of figure |3] where the total 
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Figure 4: Comparison of the lateral dependences with full bandwidth amplitudes for a 
vertical air shower with a primary energy of 10^^ eV predicted by REAS3 (left) and MGMR 
(right). The figures display the maximum absolute field strength at a given lateral distance 
of an observer to the shower axis. 



spectral electric field strength is shown as a function of frequency. For ob- 
servers with lateral distances larger than 400 m, the results match accurately 
except for the incoherent noise at high frequencies. At 100 m observer dis- 
tance, REAS3 predicts a flatter frequency spectrum than MGMR, especially 
in the frequency range of 30 — 80 MHz, in which most of the experiments mea- 
sure. The predicted amplitude of the radio emission in different azimuthal 
directions of the air shower can be compared with the lateral distributions of 
the peak amplitude of unlimited bandwidth pulses, i.e. the unfiltered signal, 
as illustrated in figure HI The lateral distributions derived from REAS3 and 
MGMR show an evident east-west-asymmetry. This asymmetry is explained 
by the time-varying charge excess in air showers causing a radiation contri- 
bution with a radial polarisation signature (cf. 0, Isij). Therefore, the two 
models give evidence that the radio emission from cosmic ray air showers is 
not purely v x B polarised. 

Looking at the contour plots of the 60 MHz spectral emission component 
displayed in figure O the same signature is visible. For radiation due to 
pure geomagnetic emission, no contribution for the north-south polarisation 
would be expected for the vertical air shower, as well as a symmetric pattern 
in the east- west polarisation component. The comparison of the contour 
plots, however, shows that there are still deviations between MGMR and 
REAS3, in particular in the strength of the east-west asymmetry predicted 
by the models. 
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Figure 5: Contour plots of the 60 MHz field strength for the emission from a 10^^ eV 
vertical air shower. From left to right: total field strength, north-south and east-west 
polarisation component. Contour levels are 0.1 //Vm^^MHz^^ apart. The closest position 
of the simulated observers to the shower core is 50 m. Upper row: REAS3. Lower row: 
MGMR 



3.2. Emission from an inclined air shower 

In this section, the resuhs for an inchned air shower with primary energy 
of 10^^ eV and a zenith angle of 50° are compared. The azimuth angle of 45° 
denotes that the shower is coming from south-east (i.e., pointing to north- 
west). Figure |6] shows the raw pulses simulated with both models for this 
geometry. Please note that the zero time corresponds to the time when the 
primary particle would hit the ground. Since the observers are located north 
of the shower core and the air shower is coming from south-east, the emission 
arrives later than in the vertical case. Close to the shower core, the predic- 
tions of REAS3 and MGMR differ nearly by a factor of three, whereas the 
results are almost the same for larger lateral distances. The pulses derived by 
REAS3 close to the shower core exhibit higher amplitudes than the pulses 
obtained with MGMR. The larger deviations close to the shower core are 
also evident in figure [71 where the lateral distributions are shown. The am- 
plitudes predicted by REAS3 increase significantly with smaller distances to 
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2000 



Figure 6: Comparison of the east-west polarisation component emitted by a 50° inclined 
air shower with a primary energy of 10^^ eV for REAS3 (left) and MGMR (right). The 
figures show pulses for observers at different lateral distances to the shower core. With 
increasing distance, the results converge. For small distances, the predictions of both 
models differ by a factor of three. 




Figure 7: Comparison of the lateral dependences with full bandwidth amplitudes for an 
air shower with zenith angle of 50° and primary energy of 10^^ cV predicted by REAS3 
(left) and MGMR (right). The figures display the maximum absolute field strength at a 
given lateral distance to the shower core for observers along various azimuthal directions. 



the shower axis. For MGMR, the lateral distribution is somewhat flattening 
to the center. Furthermore, the asymmetries between each azimuthal ob- 
server direction, i.e. north, east, south and west, are larger for MGMR than 
for REAS3. In REAS3, the observers in the north and east receive nearly the 
same signal as the observers in the south and west. For MGMR this is not 
true which might be related to the moving dipole included in MGMR but 
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Figure 8: Fraction of charge excess with respect to the total number of electrons and 
positrons as a function of atmospheric depth. With larger atmospheric depth the fraction 
rises. The large fluctuations below 200 are due to small particle statistics. 



not in REAS3. The reason for the discrepancies in the azimuthal symmetry 
might be the different implemenations of the charge excess in MGMR and 
REAS3. In MGMR, the fraction of the charge excess with respect to the to- 
tal number of electrons and positrons is assumed to be constant, whereas the 
longitudinal profile taken into account in REAS3 is changing from roughly 
14% to 26% as shown in figure M 

Above all, the differing predictions in the amplitudes might be a hint 
that close to the shower core, the details of the air shower model, which dif- 
fers in REAS3 and MGMR become important. The air shower model has a 
larger impact on inclined air showers than on nearly vertical ones, as iden- 
tical ground distances correspond to smaller effective lateral axis distances. 
Moreover, the geometrical distance between observer and shower maximum 
increases with larger zenith angle. Hence, we discuss the influence of the air 
shower model in greater depth in section HI 

3. 3. Specific magnetic field configurations 

In addition to the realistic air shower geometries shown in the previ- 
ous sections, it is interesting to look at more contrived situations such as 
some special magnetic field configurations, since the Earth's magnetic field 
is responsible for the geomagnetic radio emission in air showers. For this 
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Figure 9: Comparison of the pulses emitted by a vertical 10^^ eV air shower in the absence 
of a magnetic field (thin lines) and a magnetic field parallel to the shower axis (thick lines). 
Left: REAS3. Right: MGMR. 



comparison, the magnetic field was once switched off completely and once 
was chosen to be parallel to the air shower axis. With these special con- 
figurations of the magnetic field, the influence of the radiation due to the 
time- variation of the net charge exces^ is studied. Figure M shows that the 
results for these two configurations are indeed similar. REAS3 and MGMR 
predict both nearly unipolar pulses. This is not surprising since the longi- 
tudinal development of the charge excess component as displayed in figure 
[To] reaches the shower maximum only short before the observing height of 
1400 m which corresponds to 875 for a vertical air shower. The predic- 
tions of REAS3 and MGMR for the pure charge excess case agree very well. 
The differences in the strength of the pulses are in the order of 10-20%. The 
two models show also the expected radial polarisation pattern as displayed 
in the contour plots of figure [TTl The agreement between both models for 
the two situations shown in this section are direct evidence for the increased 
understanding of the radio emission mechanism from extensive air showers. 

4. Comparison of MGMR with simplified REAS3 

Summing up the results of the previous sections, the comparison between 
REAS3 and MGMR shows an overall agreement within a factor of ~ 2 — 3. 



^This means "Askaryan radiation" [26|,|27| but without Cherenkov-like effects since the 
refractive index is unity. 
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Figure 10: Longitudinal development of the charge excess component for the shower with- 
out magnetic field (solid red) and for a magnetic field parallel to the shower axis (dashed 
blue). 



Remembering that a few years ago the models predicted even quahtatively 
different pulse shapes and consequently differing characteristics in the fre- 
quency spectra, this agreement is a milestone in the understanding of radio 
emission from cosmic ray air showers. It should be stressed once more that 
the models are technically and conceptually very different and completely 
independent from each other. 

Looking at the details of the results from the comparison between both 
models, it is obvious that there are still deviations which are too large to 
be ignored. Especially for the inclined shower studied in section 13.21 the 
assumption that the shower model is important was made. In several figures, 
it was shown that the pulse amplitudes predicted by REAS3 are larger than 
the amplitudes from MGMR, at least for small observer distances to the 
shower core. The reason for this deviation has to be studied in detail, since 
also these discrepancies need to be understood. For the future, the aim is to 
understand the existing differences not within a factor of 2-3, but within less 
than 10%. Thus, the influence of the differences in the underlying air shower 
models has to be studied and it has to be asserted if these are responsible 
for the deviations. One possibility to study the influence of the air shower 
models is to simplify the air shower model implemented in REAS3 to get a 
more similar model to the one used in MGMR and to study the effects on 
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Figure 11: Contour plots of the 60 MHz field strength for the emission from a lO'^^eV 
air shower where the magnetic field is switched off completely. From left to right: to- 
tal field strength, north-south and east-west polarisation component. Contour levels are 
0.1 /LtVm~^MHz^^ apart. The closest position of the simulated observers to the shower 
core is 50 m. Upper row: REAS3. Lower row: MGMR 



the radio emission. In the following section, this will be done. 

4.1. The underlying air shower models 

To adapt REAS3 for comparing the details with MGMR, first the differ- 
ences of the air shower models need to be identified. One major difference is 
that the spatial distribution of the particles in the shower pancake is differ- 
ent. To mimic the arrival time distribution of MGMR in REAS3, we replaced 
the longitudinal displacement of the particles in the shower pancake with the 
F— probability distribution function given in section [TTTl 

Next, the lateral displacement of the particles in the shower pancake had 
to be changed. In MGMR, this displacement is not directly considered, but 
a radiation contribution from a static dipole with a length of 15.0 m is added 
to the overall radio signal. Hence, in REAS3, we shifted the electrons by 
7.5 m in the eastern direction of the shower axis and the positrons by 7.5 m 
in the western direction and switched off the lateral distribution of the par- 
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tides in the shower pancake. Further parametrisations used in MGMR have 
not been considered for this comparison since they influence the predictions 
less than the modification discussed here. However, they have been briefly 
discussed at the end of section [2j The results when replacing the arrival time 
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Figure 12: A comparison between the radio pulses predicted by MGMR and a modified 
version of REAS3 where the longitudinal displacement of the particles in the shower pan- 
cake follows the distribution used in MGMR and the lateral distribution is replaced by a 
systematic offset of electrons and positrons. The pulses for a 10^^ eV vertical air shower 
are shown for observers 100 m (top) and 800 m (bottom) north of the shower core. 

distribution and the lateral spread of the particles are displayed in figure [I2] 
for observers at 100 m north from the shower core and 800 m north from the 
shower core of the vertical air shower with primary energy of 10^^ eV. It is 
obvious that the modifications influence the results much stronger for small 
observer distances than for larger observer distances. At 100 m distance, the 
pulse amplitude predicted by the modified version of REAS3 deviates by less 
than 30% from the pulse amplitude predicted by MGMR whereas the original 
REAS3 pulse amplitude differed by a factor of 1.7. At 800 m, the modified 
version of REAS3 obtains the same pulse as MGMR and the original version 
of REAS3. This underlines the assumption that the details of the air shower 
model are of little significance for large lateral observer distances. This is 
reasonable since the time-scales determined by the spatial distributions of 
the individual particles in the shower pancake dominate at small lateral dis- 
tances whereas the geometrical time delays resulting from the propagation of 
the radio emission from the particles in the shower pancake to the observer 
dominate at larger lateral distances. 

As shown in section 13. 2^ the predictions of MGMR and REAS3 differed 
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Figure 13: A comparison between the radio pulses predicted by MGMR and a modified 
version of REAS3 where the longitudinal displacement of the particles in the shower pan- 
cake follows the distribution used in MGMR and the lateral distribution is replaced by a 
systematic offset of electrons and positrons. The pulses for a 10^^ eV air shower with a 
zenith angle of 50° are shown for observers 100 m (top) and 800 m (bottom) north of the 
shower core. 



the most for inclined air showers. We also used the simplified REAS3 to 
study the infiuence of the air shower model for the inclined air shower com- 
ing from south-east. The results of this comparison are shown in figure [131 
With the original versions of MGMR and REAS3, the predictions differed 
by a factor of ~ 3 for the observer 100 m north of the shower core. With the 
simplification of REAS3, the differences became much smaller, i.e. roughly 
50%. Moreover, the arrival times of the pulses were different for REAS3 
and MGMR. For the modified version of REAS3 and MGMR, they become 
similar. As expected, the infiuence of the details of the air shower model is 
less at the observer distance of 800 m but still visible and it leads to a better 
agreement between the modified version of REAS3 and MGMR. 

In general, the results of this detailed comparison showed that the differ- 
ences between REAS3 and MGMR are mostly determined by the different 
treatment of the distribution of particles in the shower front which is sim- 
plified in MGMR by a constant pancake thickness. Since we modified the 
spatial distribution of particles in the simplified REAS3 version which led 
to better agreement with MGMR, this gives a strong indication that the 
pancake thickness in MGMR is overestimated. Hence, in the following we 
investigated the differences more closely and discuss a proper value for the 
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pancake thickness used in MGMR. 

4-2. Discussion of a realistic pancake thickness in MGMR 

As stated in the introduction, the pulse shape is a direct reflection of the 
important length scales involved in an air shower. In MGMR, the pulse is 
shown to depend directly on characteristic functions describing the shower 
evolution, where the most important ones are the shower profile, the drift 
velocity, the charge-excess in the shower and the particle distributions in the 
shower front (pancake thickness). The shower profile in MGMR can either 




Figure 14: MGMR simulations for different values of the pancake thickness parameter 
L. The simulations are done for a 10"'^^ eV vertical shower with the observer positioned 
100 m to the North of the impact point. The most realistic value for the shower pancake 
in MGMR for the air showers considered in this paper is L = 2 ni (cf . fig [T2| . 



be given as a direct input from Monte-Carlo simulations, or, if the previous 
is not the parametrised profile is used. The drift velocity and the 

fraction of charge excess act as scale parameters for the geomagnetic and 
the charge excess radiation, respectively. The pancake thickness and the 
shower profile functions determine the coherence length and thus the fre- 
quency maximum of the pulse power spectrum. At large impact parameters 
the determining length scale is given by the shower profile and the pancake 
thickness is less important, while at small impact parameters this is reversed 
and the pancake thickness becomes the important parameter. With this in 



19 



mind, we can now also understand why the differences between MGMR and 
REAS3 are more prominent for small impact parameters. Both simulations 
use the same shower profile, and hence predict similar pulses for large im- 
pact parameters. As shown in this chapter, for small impact parameters the 
spatial particle distribution in the shower front becomes the important scale 
parameter. Since REAS3 simulations contain more high frequency compo- 
nents than MGMR in this regime, the logical conclusion would be that the 
pancake thickness parameter L (Eq. [3]) in MGMR is overestimated. 

In the earlier MGMR calculations jij, the pancake thickness was taken 
from the measured particle distribution in the shower front at ground, L = 
10 m, which is biased towards relatively large distances from the shower core. 



In the later version [10| used for this comparison, the average pancake thick- 
ness was obtained from Monte-Carlo simulations which resulted in L = 3.9 m. 
Within this comparison between REAS3 and MGMR and with Monte-Carlo 



simulations [28[ , we have shown that this yields an overestimate since in gen- 



eral the most important part of the electromagnetic pulse is emitted from 



distances close to the shower axis [29| where the pancake thickness is ex- 



tremely small. In the MGMR implementation for this paper the lateral 
distribution of the shower particles is ignored, and the effects are taken into 
account by an effective pancake thickness. Therefore, the value of L = 0.5 m 



as obtained in [28|, where the lateral distribution is taken into account, is 
too small and L = 2m is a more realistic value. This also follows from fig- 
ure [m where we show the influence of the pancake thickness parameter L 
on the simulated pulse shape. This is done for the 10^^ eV vertical shower 
for an observer position 100 m north of the shower core and thus has to be 
compared with figure O Since the relative importance of the real pancake 
distribution and the lateral distribution depends on the observer geometry 
the effect is estimated only on average. 



5. Conclusions 

In the present article, we discussed a detailed comparison between REAS3 
and MGMR. To minimize the influence of differences in the modelling of the 
atmosphere, this comparison has been carried out for a refractive index of 
unity. As first result, we achieved an overall agreement within a factor ~ 2 — 3 
at short distance, improving with increasing distance from the shower core. 
This agreement is already a breakthrough in the theoretical understanding of 
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the emission processes of radio emission from cosmic ray air showers because 
a few years ago, the two models predicted even different pulse shapes and 
with this deviating frequency spectra. Once again, we want to emphasise 
that REAS3 and MGMR are technically very different and completely inde- 
pendent from each other. The importance of modelling the air shower and 
with this the shower front, especially for small observer distances came up 
when discussing the inclined air shower. 

To study the influence of the air shower model on the predicted pulse 
shapes, the air shower model implemented in REAS3 was modified to get a 
more similar model to the one used in MGMR and to study the effects on 
the radio emission. We could verify that the differences in the air shower 
models are responsible for the observed deviations. 

In fact, with the modified version of REAS3 it was possible to reproduce 
the results derived with MGMR on a 20-50% level, depending on the shower 
geometry and the observer position. This indicates that the parametrisation 
of the thickness of the shower pancake used in MGMR needs to be reconsid- 
ered. A more realistic thickness of the pancake in MGMR is achieved with a 
value of L = 2 m. As this parameter L depends not only on the longitudinal 
but also on the lateral distribution of the particles in the shower pancake 
which is not explicitly taken into account in MGMR, it is difficult to make a 
hard estimate for its value. Furthermore, we showed that for larger observer 
distances, the details of the underlying air shower models are not very impor- 
tant. Hence in the first part of this comparison, the predictions of MGMR 
and REAS3 for observers far away from the shower core were already similar. 

To conclude, the comparison presented in this article showed the enor- 
mous increase of understanding the theory of radio emission from air showers 
achieved in the recent past. 
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